PREPRINT-V2-08-2012 



Converting Classical Theories to Quantum Theories by 
Solutions of the Hamilton-Jacobi Equation 



Zhi-Qiang Guo and Ivan Schmidt 

Departamento de Fisica y Centro Cientifico Tecnologico de Valparaiso, 
Universidad Tecnica Federico Santa Maria, 
Casiila 110-V, Valparaiso, Chile 

E-mail: zhiqiang.guo@usm.cl, ivan.schmidt@usm.cl 

Abstract: By employing special solutions of the Hamilton-Jacobi equation and tools from lattice theories, we 
suggest an approach to convert classical theories to quantum theories for mechanics and field theories. Some 
nontrivial results are obtained for a gauge field and a fermion field. For a topologically massive gauge theory, 
we can obtain a first order Lagrangian with mass term. For the fermion field, in order to make our approach 
feasible, we supplement the conventional Lagrangian with a surface term. This surface term can also produce 
the massive term for the fermion. 

Keywords: Dc Dondcr-Wcyl Theory, Topologically Massive Gauge Theory, Mass Generation for Fermion 



Contents 



1 Introduction 2 

2 Examples of Mechanics 2 

2.1 Linear Theories: Harmonic Oscillator 2 

2.1.1 Path Integral Quantization for the Harmonic Oscillator 2 

2.1.2 Hamilton- Jacobi Equation for the Harmonic Oscillator 3 

2.1.3 The Pauli Short-Time Kernel 4 

2.1.4 Deriving Discretized Lagrangians from Special Solutions of the Hamilton- Jacobi Equation 5 

2.2 Nonlinear Theories: Double- Well Potential 6 

2.2.1 Tools for Solving the Nonlinear Hamilton- Jacobi Equation 6 

2.2.2 Derivation of the Lattice Lagrangian 8 

3 Examples of Scalar Field Theories 9 

3.1 Linear Theories: Massive Scalar Theories 9 

3.1.1 Covariant Hamilton- Jacobi Equations 9 

3.1.2 Solutions of Covariant Hamilton- Jacobi Equations 10 

3.1.3 Derivation of the Discretized Lagrangian 11 

3.2 Nonlinear Theories: Scalar Theories with \<p A Potential 13 

3.2.1 De Donder-Weyl Equation and Its Solution 13 

3.2.2 Derivation of the Discretized Lagrangian 15 

4 Topologically Massive Gauge Theory 15 

4.1 De Donder-Weyl Equation for Gauge Theories 16 

4.2 Derivation of the Discretized Lagrangian 18 

5 Mass Generating Mechanism for a Fermion Field 20 

5.1 De Donder-Weyl Theory for a Fermion Field 20 

5.2 Derivation of the Lattice Lagrangian for a Fermion field 23 

6 SU(2) Yang-Mills Gauge Theory 23 

6.1 De Donder-Weyl Equation of the Yang-Mills Gauge Theory 23 

6.2 Derivation of the Lattice Lagrangian 26 

7 Further Discussions 27 

7.1 Discussions on Solution Dependence 27 

7.2 comparisons with conventional lattice gauge theories 28 

8 Conclusions 28 



- 1 - 



1 Introduction 



Quantum theories have achieved tremendous success in the passed century. There are two conventional ap- 
proaches to convert classical theories to quantum theories: canonical quantization and path integral quan- 
tization. These two approaches employ two basic objects in classical theories: the Hamiltonian and the 
Lagrangian, which are two equivalent tools to formulate classical theories. However, besides these, we know 
there is a third equivalent way to formulate classical theories: the Hamilton- Jacobi equation. In this paper, 
we show that it is possible to derive quantum theories from classical theories by employing special solutions of 
the Hamilton- Jacobi equation, together with some tools from lattice theories. This third approach turns out 
to be consistent with the path integral approach for most of the cases. However, we also can obtain several 
new sectors for gauge fields and fcrmion fields. Specifically, we can obtain a massive Lagrangian of the first 
order for the topological massive gauge theory introduced in [1, 2]; While for a fermion field, we find that the 
mass of a fcrmion can be produced by a surface term, which is a mass generating mechanism similar to the 
topological massive gauge theory. 

This paper is organized as follows. In section 2, starting with classical mechanics, we introduce the 
basic tools and methodology used to derive quantum mechanics from special solutions of the Hamilton- Jacobi 
equation. We also introduce tools to find solutions of the Hamilton- Jacobi equation for nonlinear theories in 
section 2. We turn to scalar field theories in section 3. In this section, we introduce the covariant Hamilton- 
Jacobi equation for field theories, that is, the De Donder-Weyl approach [3-5] for field theories. Along with 
the discussions regarding mechanics, we introduce tools for solving the Dc Donder-Weyl equation for nonlinear 
field theories, and stress the differences between mechanics and field theories. In section 4 and section 5, we 
discuss the topological massive gauge theories and fermion fields separately. Section 6 is devoted to SU(2) 
Yang-Mills theories. We only obtain restricted solutions for Yang-Mills theories, and no firm conclusions can 
be drawn from these restricted solutions. We provide further discussions and conclusions in sections 7 and 8. 

2 Examples of Mechanics 

We discuss mechanics in this section. Take the harmonic oscillator, for example. The logic structure of these 
discussions is as follows: First we present the path integral quantization for the harmonic oscillator and its 
lattice definition: Then we display its Hamilton- Jacobi equation and find several solutions for the Hamilton- 
Jacobi equation. Based on the discussions above, we can find a close relation between the lattice definition 
of path integral quantization and a special solution of the Hamilton- Jacobi equation. This relation will make 
it feasible to convert classical theories to quantum theories by the special solution of the Hamilton- Jacobi 
equation. 

2.1 Linear Theories: Harmonic Oscillator 

2.1.1 Path Integral Quantization for the Harmonic Oscillator 

For a massive harmonic oscillator in one dimension, its Lagrangian is given by 



where x — ^ is defined. For a time evolution from t a to tb, the Green function or the Feynman kernel is 
defined by the path integral in configuration space 



where AT is a normalization factor. Here and hereafter we use natural units, so the Planck constant H = 1 is 
assumed. Discretizing the interval [t a , tt,] into N equivalent smaller intervals with length e = tb jv *° , that is, 



L(x,x) = - 



•2 2 2 

mx muj x 

2 



(2.1) 
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let t a = to, t\ — to + e, ■ ■ ■ , tk = to + fee, • • • , t b = tjsi, then the path integral (2.2) can be regarded as the 
limit of multi-integrals and can be expressed as the following lattice version 



AT-lp , . i , N-l 

,2( x j+l ' X 3 



K(xb,t b ;x a ,t a ) = lim / dx 3 ( — — ) 2 explie V* — ( -^±i -i 

N^toa J -- L / Kline/ z — ' 2 LV e 

e^O J=l L,/ k j=l 

For infinitesimal time evolution e, we have the short-time Feynman kernel 



— w 



(2.2) 



^W, + e;^-) = (^ (2.3) 

So the Feynman kernel (2.2) of finite time evolution also can be regarded as multi-convolutions of short- 
time Feynman kernels (2.3) 

K(x b , t b ; x a , t a ) = Yira J I ■ ■■ J dx N ^ ■ ■ ■ dx 2 dx 1 

K(x b ,t N -! + e;x N -!,t N -i) ■ ■ ■ K(x 2 ,h + e;x 1 ,t 1 )K(xi,t a + e\x a ,t a ). (2.4) 

The Gaussian integrals in (2.2) can be performed in sequence and we can get a closed formulation for the 
Feynman kernel 

( m0J \ 5 imuj . o 2 1 1 

K{x b ,t b ;x a ,t a ) = I — — : 7- —r) exp< — —[cosu}{t b -t a )(x b +x a )-2x b Xa\>. (2.5) 

\2nism.oj{t b — t a )J y 2smui{t b — t a ) J 

2.1.2 Hamilton-Jacobi Equation for the Harmonic Oscillator 

The Hamilton-Jacobi equation was independently introduced by Hamilton and Jacobi from different ap- 
proaches. In this section, we only give a pedagogical introduction based on the independent integral of 
Hilbert [6], which has been applied to field theories by De Donder [3] and Wcyl [4, 5]. We caution that our 
introductions only work well for regular Lagrangians, which we always work with in this paper. For non-regular 
Lagrangians, we refer the reader to the more rigorous discussions in [4, 5, 7-9]. 

For a system of mechanics in one dimension, its classical aspects can be formulated by a Lagrangian 
L(q(t),q(t),tJ, where q(t) is the coordinate and q = ^| is defined. After performing the Legendre transfor- 
mation p — we can get the corresponding Hamiltonian H(p, q, t) = qp — L. The Hamilton-Jacobi equation 
can be derived as follows. Suppose there is a function S(q, t) that depends only the coordinate but not its time 
derivative q(t), then the Lagrangian can be regarded as an independent integral of Hilbert with the following 
meaning 



which also can be expressed as 



dS _ OS dq dS 



Designating p = ^ , we can get the Hamilton-Jacobi equation straightforwardly 

dS dqdS T dS dq T dS TT ( dS \ n ,„ n , 

m + iTq- L =m + i p - L =m +H { q >8-q>V= - (2 ' 8) 

For the harmonic oscillator (2.1), its Hamiltonian is given by 

H = J-p 2 + \ra 2 uj 2 x 2 . (2.9) 
2m 2 
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Designating p = we derive its Hamilton- Jacobi from eq. (2.8) as 



dS 1 /dS\ 2 1 2 2 „ , B<n , 

— + — — +-mw¥ = 0. 2.10 

dt 2m\dxJ 2 y ' 

We can find solutions of eq. (2.10) following two different approaches as follows: 
Type-(I): Assuming the solution is a polynomial of x, we have 

S= \f{t)x 2 + h(t.)x + g{t). (2.11) 

Substituting this assumption into eq. (2.10), and letting the coefficients of x to be zeros, we can get equations 



ldf(t) , 1 . u . 3 , 1 2 



2 ■^/(') J +2 mW =°' (2 ' 12) 

« + I,(t)*(t) = 0, (2.13) 

f + ^#(t) =0. (2.14) 

Solving these ordinary differential equations (ODEs), we can get the solution for S to be 

(la): £-—[cosLu(t-t a )(x 2 +x 2 )-2xx }, (2.15) 

2 sinu;(r — to) 

(lb): Si = ^ [-smLj(t-to)(x 2 +x 2 )-2xx Q ], (2.16) 

Z COS OJyt — to J 

where xo, to are integration constants. 
Type- (II): Assuming that 

S = -E(t-t ) + W(x), (2.17) 
we can get a solution which we are familiar with in classical mechanics 

S = -E(t -t ) + — arctan ( . mU}X ) + \x\j2mE - m 2 oj 2 x 2 , (2.18) 

V 7 m \V2mE-m 2 Lo 2 x 2 J 2 1 ' 

where E, to are constants of integral. 
2.1.3 The Pauli Short-Time Kernel 

Based on the discussions above, we should have noticed the resemblance between the Feynman kernel (2.5) 
and the solution (2.15) of the Hamilton- Jacobi equation. Of course, this resemblance is not so surprising. The 
reasons are as follows: For a Lagrangian of quadratic interaction, its Feynman kernel can be given exactly by the 
Wcntzcl-Kramcrs-Brillouin semiclassical approximation; While the Wcntzcl-Kramers-Brillouin approximation 
involves the classical action of the system. Therefore, the resemblance between (2.5) and (2.15) only holds 
for quadratic interaction or for linear theories. For nonlinear theories, there would be no such resemblance. 
However, we will show that a connection between the discretized Lagrangian and special solutions of the 
Hamilton- Jacobi equation can still be constructed, not only for linear theories but also for nonlinear theories. 
In quantum mechanics, the link connecting them is Pauli's formula or Pauli short-time kernel [10, 11]. 
For infinitesimal time evolution from t to t + e, the Pauli short-time kernel is defined by 

K,(,.. + ^..)-(i)*(-^) 1 ^>. MO 

Pauli proved that this kernel will satisfy the same equation as the Feynman short-time kernel (2.3) if the 
function S(q,q';e) satisfies the Hamilton- Jacobi equation 

+*(,,§?)_ a (2.20) 
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For the harmonic oscillator, the solutions of its Hamilton- Jacobi equation are given by eqs. (2.15), (2.16) and 
(2.18). Obviously the Type-(Ia) solution eq. (2.15) is suitable to formulate the Pauli short-time kernel while 
the solutions Type-(Ib) (2.16) and Type-(II) (2.18) arc not. The Pauli short-time kernel for the harmonic 
oscillator is therefore given by 

/ moj \ 2 f m ui r 2 2\ 1 1 

Kp[x j+1 ,t s+ i;x j ,t j ) = ^ — ainu} J ex PVYSn^r° Sa;e ^ +1 +Xj ' ~ 2 ^' +1 ^J f {2 ' 2l > 

Here we use the subscript P to differentiate the Pauli short-time kernel from the Feynman short-time kernel. 
The convolution of the Pauli short-time kernel (2.21) has the semi-group property, that is, 

Kp(x j+2 ,tj + 2e;Xj,tj) = J dx J+ iK P (x j+2 , t j+ i + e; xj+i, tj +1 )K P (xj +1 , tj +e;Xj,tj). (2.22) 

So its multi-convolutions will give the same results as the Feynman kernel (2.5). The calculations of the multi- 
convolutions of the Pauli short-time kernel (2.21) are simple and straightforward; While the multi-convolutions 
of the Feynman short-time kernel (2.3) are very complicated. 

2.1.4 Deriving Discretized Lagrangians from Special Solutions of the Hamilton- Jacobi Equation 

The discussions in section 2.1.3 suggest us an approach to convert classical theories to quantum theories. The 
logic is as follows: Employing the multi-convolutions of the Pauli short-time kernel (2.21), we can get the 
same results as that we get with the Feynman short-time kernel (2.3); While the Pauli short-time kernel (2.21) 
can be determined by the special Type- (la) solution eq. (2.15) of the Hamilton- Jacobi equation. The logical 
deductions above can be reversed: we can begin with the Type- (la) solution eq. (2.15), then we can get a 
discretized version of the path integral via the multi-convolutions of Pauli's formula (2.19). 

The equivalence of the results of these two approaches can be understood in another way. Expanding the 
Pauli short-time kernel (2.21) with the small parameters e and keeping the terms of leading order, we get 

2 [cos uie(x 2 +1 + x 2 ) — 2xj + \Xj\ = (1 + cosLue)(xj +1 — Xj) 2 — (1 — cosa;e)(a;j + i + Xj) 2 (2.23) 

(2.24) 







< m n 


\2iri smujeJ 




K2-Kit) 



m ui 
cxp< i- 



2 sin we 



cosuie(x 2 +1 + x 2 ) — 2xj + ix 3 



> cxp< le 



m /Xj+i - x 3 \ 2 1 2 ( x 3+i + x j\ 2 ] \ 



(2.25) 



What we should realize from this approximation is the discretized Lagrangian in the lattice definition (2.2) of 
the path integral is recovered from the Type-(Ia) solution eq. (2.15) of the Hamilton- Jacobi equation in the 
approximation of small e. 

The derivation of the discretized Lagrangian can be generally summarized in the following procedure. 
For a solution of the Hamilton- Jacobi equation S(x, t; xq, to) with some constants (xo,to) which need to be 
designated by the initial conditions and boundary conditions, its corresponding discretized Lagrangian can be 
defined by 

f _ S(x b ,t b ;xo,to) - S(x a ,t a ;x ,t ) 

-^lattice — " 1 ■ (Z.ZO) 

% — ta 

Here [t a , is the lattice interval. The constants (aio,£o) need to be designated appropriately. Substituting 
the Type-(Ia) solution cq. (2.15) into eq. (2.26), we get 

t 1 ( m Ld \/2 2\r.l 

^lattice = -riTT- 7T c osu;(ib - t )(x b + x Q ) ~ 2x b x 

t b - t a \ 2 smw(t(, - t ) 

- n T^[™Mta-to)(x 2 a +x 2 )-2x a x ] \. (2.27) 

Z smu}{t a — to) J 
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Taking the limit xq — > x a , then the limit to — > t a , we get 

^lattice = lim lim Lattice = 7 1 ^ | 7T ~ 7T 7T [cOSw(t b - to) (^6 + x a) ~ ^ b X a ] 

For infinitesimal lattice spacing t b —t a = e, following the procedure in eq. (2.25), eq. (2.28) can be approximated 
as 

lfm U! . . . , 2 2\ r. 1 

^lattice = 7 1 TT ~ 71 -tI cosw (*& - MO^b + X a ) - 2x b X a \ 

t b — t a I Z smuj(t b — t a ) 
tb -^ ™(*b- x «) 2 l rnu >( x » + x °)\ (2.28) 



2 V t b — t a ) 2 
The lattice Lagrangian is recovered again. 

2.2 Nonlinear Theories: Double- Well Potential 

In this section, taking the double-well potential for example, we deal with nonlinear theories. For nonlinear 
theories, the Hamilton- Jacobi equation is generally difficult to solve. We introduce the tools to handle the 
nonlinear theories and show that the discussions in section (2.1) can also apply to nonlinear theories. 

2.2.1 Tools for Solving the Nonlinear Hamilton- Jacobi Equation 

For a particle in the double-well potential, its Lagrangian is given by 

L(x, x, t) = l -mx 2 - l -g 2 (x 2 - v 2 ) 2 , (2.29) 
where v is constant. After Legendre transformation p = we get its Hamiltonian 

H{x lP ,t) = ±-p 2 + 1 -g 2 {x 2 -v 2 ) 2 . (2.30) 
The canonical Hamiltonian equations of motion are 



dH _ p 
dp m ' 
dH 1 
~dx ~ ~T 



(2.31) 

P = -^ = ~ffM* a -« a ). (2-32) 



While the corresponding Eulcr-Lagrange equation of motion is 

1 



mx + ^g 2 x (x 2 - v 2 ) = 0. (2.33) 



Designating p = , its Hamilton- Jacobi equation is 



dS 1 /<9SV 1 



^L + J^I^L) +-g 2 ( x 2 -v 2 ) 2 =0. (2.34) 
dt 2m \ dx J 8 y K 1 K ' 

This nonlinear equation is difficult to solve. However, employing the "embedding method" introduced in [12], 
we can find a series solution for eq. (2.34); While this series solution is enough to satisfy our purpose in this 
paper. 

The "embedding method" is as follows. Suppose that we seek a series solution of the type 

S(x, t) = S*(t) + A(t) [x - /(<)] + R(t) [x - f(t)} 2 

+ K(t) [x - f(t)f + M(t) [x - /(t)] 4 + N(t) [x - f{t)f + ■ ■ ■ , (2.35) 



- 6 - 



where f(t) is a function which will be given later. So it seems like we expand S(x,t) in a series around 



a function fit) in eq. (2.35). The potential function V(x) 
expanded by the following identities 



b 2 (* 2 



is a polynomial, which can be 



v(x) = i 9 V - « 2 ) 2 = \g 2 \(x fit) + f(t)f - v 2 



(x - f(t) + fit)) 2 = [x- fit)} 2 + 2f(t) [x - fit)} + f{tf 



(2.36) 



Notice that we expand these functionals around a function f(t) but not 0, so the combination [x — f(t)] always 
remains. Substituting the expansions (2.35) and (2.36) into the Hamilton- Jacobi equation (2.34), collecting 
the terms of [x — f(t)] of the same power, and letting the coefficients of this series to be zeros, we get 



x - fit) 
x - f{t) 

x - /(*)' 
x - f{t) 
x - fit) 



dt 
dA(t) 
dt ' 



A 



dfjt) 
dt 



dRjt) 
dt 



- 3H 

dK(t) 
dt 



2R(^ 

dfit) 
dt 

-AM 



dfjt) 
dt 
A 



A 
m 



1 



gifit) [fitf-J 



m 
dfjt) 
dt 



-R 2 + \g 2 [lf(t) 2 -v 2 ] 

vn 4 J 



.4 
m 



6 



-KR 



1 



9 2 fit) 



dM(t) 
dt 



-5N 



dfjt) 
dt 



mi im 



l&MR) + -g 2 
8 



0. 


(2.37) 


0. 


(2.38) 


0. 


(2.39) 


0. 


(2.40) 


0. 


(2.41) 



where we only display terms of power not higher than 4 while those of higher power are omitted for convenience. 
These ODEs are generally difficult to solve. The key ingredient of the "embedding method" is we can use 
some assumptions to simplify these ODEs. We suppose 



dm 

dt 



A 

- = 0, 

m 



^ + \g 2 fit) [fit? - v 2 ] 



o. 



(2.42) 
(2.43) 



then eq. (2.38) is satisfied. We noticed that equations (2.42) and (2.43) coincide with the canonical Hamiltonian 
equations (2.31) and (2.32) if we make the replacements x — > fit) and p — > A. So f{t) satisfies the Euler- 
Lagrange equation automatically 



1 



mf + -g 2 f{f-v 2 )=0. 



(2.44) 



Eq. (2.37) can be transformed as 

dS*(t) 
dt 



A 



dfjt) 
dt 



J-A 2 --g 2 \f(t) 2 
2m 8 y u v ' 



L*(fJ,t). 



(2.45) 



We see under the assumptions (2.42) and (2.43), the righthand side of (2.45) coincides with the Lagrangian 
evaluated on the configuration fit). While eqs. (2.39), (2.40) and (2.41) are simplified to be 



dR(t) 
dt 



dMjt) 
dt 



- —R 2 

m 

dKit) 



:9 2 [W? 



dt 



-KR 



;9 2 fit) 



1 

2m 



(9K 2 + 16MR) 



8 ; 



0. 
0. 

= 0. 



(2.46) 
(2.47) 
(2.48) 



We see that eq. (2.46) only includes coefficients of power less than 3, eq. (2.47) only includes coefficients of 
power less than 4, and eq. (2.48) only includes coefficients of power less than 5. So coefficients of higher power 
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terms decouple from that of lower power terms. This is a great advantage of the "embedding method". From 
the discussions above, we learned an important point in the "embedding method" , which is we need to know 
a solution for the original Hamiltonian canonical equations or the Euler-Lagrange equation. This solution 
of the equation of motion will help us find a series solution for the Hamilton- Jacobi equation. Obviously, 
this "embedding method" is useless if our purpose is to find the solution of Hamiltonian canonical equations. 
However, here the "embedding method" just satisfies our purpose because we try to find a solution of the 
Hamilton- Jacobi equation. 

A solution of eq. (2.44) can be expressed by the Jacobi's function 



f(t) = xoJacobiDN^A(t-io),fc), A 



9Xq 



2(l-(^ 

\x 



(2.49) 



where we have used the initial condition /(to) = xq. Substituting the expression of /(f) into eqs. (2.46), (2.47) 
and (2.48), we can obtain solutions for R(t), A(t) and M(t). However, because of the limiting procedures in 
eqs. (2.28) and (2.28), the solutions for small (t — to) are enough for our purpose. For small (t — to), /(f) can 
be replaced by x in eqs. (2.45), (2.46), (2.47) and (2.48). The solutions for S*{t), R{t), A(t) and M(t) can 
be given by 



S*(t) 
R(t) 
A(f) 
M(t) 



ml m 
2 r 6 

C*2 C2 to C\ 9 „, o, 

-T + — T + — " 2 r + 0{r% 



2 r 

9 c 2 a 



J 2 

2m r 5 



3w G? 



2m r 3 



2uj C 3 

Y^ 2 



24m 



-9d + 7C 2 )C 2 - 
r 



1 
720 



a; 
m 



153C*i - 56C 2 C 2 + 18g 



r + 0(r 2 ), 



2w J 



-C 3 



(2.50) 
(2.51) 
(2.52) 
(2.53) 



where r = t-t ,u= 2 [3xg - v 2 ], C 2 = Ci 
derived from eq. (2.42) straightforwardly 



A(t) 



1 



^zXq and C3 is a constant. The expression for A(t) can be 

(2.54) 



lg 2 ( x 2 - v 2 )r + 0(r 2 ) 



2.2.2 Derivation of the Lattice Lagrangian 

We have obtained a series solution for the Hamilton- Jacobi equation in section 2.2.1. Then following the 
limiting procedures in 2.1.4, we can derive the lattice Lagrangian. Because the additive properties of the 
limiting, we can implement the limiting procedures term by term. The results are given by 



S*( Xb ,t b ) S*(x a ,t a ) = _1 2 / 2 2)2 (255) 

«o-**. t b -t a g y \ a 

i— [A(t h ) (xj, - /(f fc )) - A(t ) (*„ - /(t a ))] = -\g 2 {x 2 a - v 2 )(x b - x a ),(2.56) 



lim 



lim 

tb - ta 



R{tb) {X b - /(ffc)) 2 - R{t a ) (Xa - /(fa))' 



lim 

1 



A' (t b ) (x 6 - /(t & )) 3 - A(t a ) (z a - /(f a )) 3 J = - — 9 2 a; a (a; 6 - x a ) 3 ,(2.58) 



lim 

*«• - *« 



M(t b ) (.T fc - /(tfc)) 4 - M(f a ) (z a - /(t a )) 4 



40 



g 2 (x b - x a ) 4 (2.59) 
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Here we caution that we take the limit x$^x a at first, then take the limit to—>t a . In the limiting procedures 

2 

above, we have set the constants of integral C2 and C3 to be zeros, so C\ = — -^jxo, in order that the limits 
can be well defined; Otherwise the limits will be singular function of (to — t a ). The lattice Lagrangian is then 
found to be 



^lattice 



v -- {Xa -v){x h -x a ) + - (2.60) 



2 



5 2 3 5 2 4 

Xa(,Xb X a ) 40 ~t~ ' ' ' • 

Again for infinitesimal lattice spacing % — t a = e, cq. (2.60) can be further approximated as 



j tb-rtg m (x b - x a ) 2 g 2 

-^lattice ^ „ /, , s 2 o 

x b ^x a I [t b — t a ) z O 



Xb + .X a \ 2 o 



1 2 



(2.61) 



which is the expected lattice version for the double- well Lagrangian (2.29). We saw that the terms of power 
higher than 3 disappear in the limit x — >■ x a . For the double- well potential, we can verify this phenomenon 
term by term. Thereafter, when we find series solutions for nonlinear theories, we always suppose that this 
phenomenon holds. It seems to be plausible despite being difficult to prove. 

3 Examples of Scalar Field Theories 

In this section, we discuss scalar field theories. We try to apply the discussions of section 2 to field theories. 
To make these applications feasible, some new tools are required. These new tools are the covariant Hamilton 
theories or the Dc Donder-Weyl theories for field theories, based on the multi-parameter generalization of 
Hilbert's independent integral. Taking scalar field theories for example, we will give a pedagogical introduc- 
tions of the De Donder-Weyl theory. For more details, see [4-8]. Based on these covariant Hamilton- Jacobi 
equations, we suggest an approach to derive discretized Lagrangians for field theories, similar to the procedures 
in section 2. 

3.1 Linear Theories: Massive Scalar Theories 
3.1.1 Covariant Hamilton- Jacobi Equations 

For a massive free scalar field in four dimensional Minkowski space-time, its Lagrangian is given by 

& = lird^(x)d„<f>(x) - ^m 2 ct>{x)ct>{x). (3.1) 

Here and hereafter we use the Lorentzian metric 77 = diag(l, —1,-1,-1). Following Dc Dondcr and Weyl, the 
Legendre transformation is performed with the manifest Lorentz covariance 

^ = dW¥) = d ^ {xl (3 - 2) 

which is different from the Legendre transformation only performed in the time derivative of a field; While 
tt^ can be regarded as the covariant conjugate momentum. Then the corresponding covariant Hamiltonian is 
found to be 

J? = ir^<t){x) - jSf = -t^tt" + -m a <t>(x)<t>(x). (3.3) 
From this covariant Hamiltonian, we can get the Hamiltonian canonical equations 



d»<Kx) = ^ = tt", (3.4) 
3/% = _£^L = _ TO 2^. (3.5) 
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Obviously these canonical equations imply the following Euler-Lagrangian equation 



d l *d fl (j) + m 2 <j) = 0, (3.6) 

as we expected. 

The covariant Hamilton- Jacobi equation or the De Donder-Weyl equation can de derived as follows. Sup- 
posing the Lagrangian Jz? is an independent integral of Hilbcrt, that is, JC can be expressed as the total 
derivative of a vector S 7 ' 

, _ dS^ dS^ — , , \ dS^ /„ _. 



which also equals to 



dxf * v 'defy 

Designating 



d^(x)^--J? = Q. (3.8) 



we then get the covariant Hamilton- Jacobi equation 

— + ^(W = — )=0. (3.10) 

As in classical mechanics, if we can find a complete solution for from the Hamilton- Jacobi equation (3.10), 
then we can get a solution for the Euler-Lagrange equation (3.6). However, there are some differences between 
field theories and classical mechanics. The reasons are that a classical mechanical system only depends on a 
single evolution parameter, the temporal variable t; While a system of fields depend on the variables of time 
and space dimensions. Therefore, there are some integrability conditions involved for field theories. These 
discussions can be understood as follows. If we find a solution for S^, then we can determine by eq. (3.9); 
While on the other hand, is a conjugate momentum determined by the Legendre transformation (3.2), so 
it can be expressed as the total derivative of fields. Therefore tt^ determined by eq. (3.9) is connected to the 
total derivative of fields. The foregoing are the reasons why the integrability conditions appear. For the free 
massive scalar field in this section, the integrability conditions are straightforwardly found to be 



diTfj, _ d-K v 
dx v dx^ 



(3.11) 



or more specifically 



+ ^flf =dxT + ^ {X) W (3 - 12) 



Here we noticed that tt^ is determined by eq. (3.9) but not eq. (3.2), that is, 7r M is derived from the functional 
S^. So these integrability conditions are basically some restriction conditions on S^. For more detailed 
discussions on integrability conditions, see [5, 8, 12]. 

3.1.2 Solutions of Covariant Hamilton-Jacobi Equations 

For the Lagrangian (3.1) of a massive free scalar field, its De Donder-Weyl equation is given by 

dS» 1 dS" dS u 1 2 , . , , . , . „ , 
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We employ the "embedding method" introduced in section 2.2.1 to find its solution. In order to do that, we 
need a solution of the Euler-Lagrange equation (3.6). We can find two types of solutions of equation (3.6): 



I 



(771 \ 
_rj, A = fc M fc M , r = k M (x - zY, 



II: ip(x) = (p{z) =BeecslJ(l, m-v/f), £ = (x — z)Jx — z) A 



(3.14) 
(3.15) 



where fc M and z M are constant vectors and we have used the initial conditions to normalize the solutions. These 
two types of solutions can both be used in the "embedding method" . We take the type (I) solution for example. 
Because eq. (3.13) only contains quadratical terms, we can suppose that is given by 



S" = /"(x) 4>(x)-ip{x) + h^(x) (j>(x)-(p(x) + S*»(x). 



(3.16) 



Substituting this expression into the Dc Dondcr-Weyl equation (3.13), we can get a polynomial of [<fi(x) — ip(x)]. 
Letting the coefficients of this polynomial to be zero, we get 



<t>{x) - ip(x) 
(x) - ip(x) 
(x) - <p(x) 

Similar to that in section 2.2.1, we can suppose 



d^ST" - dtftohf + + ^m 2 <p 2 (x) = 0, 



8^ + 2f»[-d ll <p(x) + / V ] + m*<p(x) = °> 

1 



» A- -7 



0. 



dh>* 



d^x) 



= -™V(«), 



(3.17) 
(3.18) 
(3.19) 

(3.20) 
(3.21) 



then eq. (3.18) gets satisfied. Obviously, Eqs. (3.20) and (3.21) are just the Hamiltonian canonical equations 
(3.4) and (3.5). They are consistent because <p(x) is a solution of the Euler-Lagrange equation. While eq. (3.17) 
is transformed to be 



dS*>* 
dxf 1 



(3.22) 



whose right-hand side is the Lagrangian evaluated at ip(x). The solutions for and S* 1 " 1 can be given by 



1 m ( m 
-— =cot — =r 

2 Va vVa ' 



Fix) 



^--sin^, 



(3.23) 
(3.24) 
(3.25) 



where r has been defined by cq. (3.14). 

3.1.3 Derivation of the Discretized Lagrangian 

In order to derive the discretized Lagrangian for field theories, we need a formula similar to eq. (2.26). However, 
some differences emerge in field theories, because what we have in field theories is the vector functional S 1 *. So 
we need a multi-dimensional exploration for the formula (2.26). Inspired by the designation of the independent 
integral of Hilbcrt, that is, cq. (3.7), we suggest the following definition for discretized Lagrangian 



lattice 



E 

i=0 



yl _ X l 



(3.26) 
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where 

sV)-sV) 



y 1 — x 1 



y 1 — x 1 



S 1 (a; o , 2 / 1 ,a; 2 ,x 3 ,0(a; o ,y 1 ,a; 2 ,x 3 )) - S 1 (x°, x 1 , x 2 , a; 3 , </>(x° , x 1 , x 2 , x 3 )) 



(3.27) 



The definitions for other indices follow similarly. Obviously, in the limit of y % — x % — > 0, this definition is just 
the independent integral of Hilbert (3.7). For the solutions in eqs. (3.23), (3.24) and (3.25), employing the 
definition in (3.26), then taking the limits of z % — > x % and tp(z) — > 4>{x), we can get 



latlii 



lim lim lattice 

3 



(3.28) 



1 m 



1 m 



/2m 



Va 



] k 



Here we have supposed the symmetrical lattice spacing, that is, y 1 — x l = e. Furthermore, supposing infinites- 
imal lattice spacing e — > and 4>{y l ) — <\>(x % ) — > 0, we get the final version of the discretized Lagrangian 

3 



■XI 



■ Acc = lim lim -lattice = Sgn(i)i — <j){y % ) 



1 J 

— m 

2 



Here we have defined the function 



sgn(i) 



i = 0, 

i = 1,2,3. 



(3.29) 



Obviously, cq. (3.29) is the lattice version of the Lagrangian (3.1) as we expected. The appearance of the 
sgn(i) is because we work on Minkowski space-time with a Lorentzian metric. For infinitesimal e — > 0, we can 
make the replacement 



4>(y l ) - 4>{x l ) e ^o 



(3.30) 



(3.31) 



> di<f>(x), 

t 

then the lattice Lagrangian will approximate to the continuous one 

-lattice i ^ ^qi iV d ll ^{x)d v <i>{x) - ]^m 2 <j)(x)<f){x)i 

which is of course the Lagrangian (3.1) we began with. 

Moreover, we should mention that the Lagrangian (3.29) is not the only lattice Lagrangian we can obtain 
by following the foregoing procedure. The reason is that we can obtain more solutions besides the solutions 
in equations (3.23), (3.24) and (3.25). Actually, we can get two more solutions of eq. (3.19) for / M 

1 m 



/f(z) 



2^ 



tanh 



VA 



r fc". 



1 m BesselY(2, m-y/^) 



{x-zY, 



(3.32) 
(3.33) 



2 BesselY(l,m-v/f) 

where r and £ have been defined in eqs. (3.14) and (3.15). Associated with solutions (3.24) and (3.25), 
eqs. (3.32) and (3.33) both construct new solutions for S 1 *. Following the limiting procedures as we just did 
above, we can derive two new Lagrangians from these two new solutions. They are given by 

1 



fS ■ = -\iTd^{x)d v 4>{x) - im 2 0(x)0(x). 



(3.34) 
(3.35) 
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The first new one has no kinetic term because /f does not contribute in the limit of e — ¥ 0; While the second 
new one has the minus kinetic term, that is, a ghost kinetic term because f£ yields the minus kinetic term. 
We might expect the integrability condition introduced in (3.11) to kill these two new solutions, but they both 
satisfy the integrability condition. We will see that this phenomenon happens for linear theories universally. 
So far we can not figure out they could correspond to two new sectors of quantum theories we can derive from 
solutions of the De Donder-Weyl equations, or they mean the procedure we have employed is incomplete so 
we need some more criteria to select the physical solution. 

3.2 Nonlinear Theories: Scalar Theories with \<fi 4 Potential 

In this section, we discuss nonlinear field theories. These discussions are in conjunction with that of linear 
theories. The differences are we can get exact solutions for linear theories, but we can only get series solutions 
for nonlinear theories. 

3.2.1 De Donder-Weyl Equation and Its Solution 

We consider a scalar field theory of Xcf> 4 potential in four dimensional Minkowski space-time, its Lagrangian is 
given by 

% = \^d^{x)d v <j>{x) - ^m^ix^x) - ±<f(x) - A, (3.36) 
where we include a density of vacuum energy A for generality. Performing the Legendre transformation 



we get the covariant Hamiltonian 



= d^(x), (3.37) 



JP = X -ix^ + im 2 0(x)0(x) + ^0 4 (x) + A. (3.38) 



The corresponding Hamiltonian canonical equations are 

dir 



&><t>(x) = -£- = tt", (3.39) 
d^^-^^~m^-^(x). (3.40) 



They imply the following Eulcr-Lagrangian equation 



^3^0 + m 2 ^ + 7^ 3 (a0 = 0. (3.41) 



In eq. (3.38), designating 



we get its De Donder-Weyl equation 



90 ' 



(3.42) 



W + lltylk + 2 m20(x)0(;E) + 4! 04(a;)+A = O - (3 - 43) 

Following the procedures in section 3.1.2, we employ the "embedding method" to find the solution of 
eq. (3.43). As in section 3.1.2, we need a solution of the Euler-Lagrange equation, which can be given by 



(p(x) = ^(z)JacobiDN^^^ r,fcj , fc = ^2 + ^^, r = p„(s - a=p^. (3.44) 
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Here p^ and are constant vectors, and we have used ip(x)\ z = <fi(z) to normalize the solution. According to 
the "embedding method" , we suppose a series solution for 



S" = S***(x) + P M (x) [<^(aj) - p(a;)] + R^(x) [4>(x) - (f(x)} 2 

+ K»{x) [<(>(x) - p(x)} 3 + M"(x) y>(x) - ^(x)} 4 + N»(x) [<f>{x) - ip{x)f 



(3.45) 



Substituting this expression into eq. (3.43), we get a series expression of [4>(x) — ip(x)]. Supposing the coeffi- 
cients of this series to be zeros term by term, we get 



[<j>(x) 


- Viz)} 


\4>(x) 




W) 


- <P(x)} 2 


\4>{x) 


- v( x )? 


[4>(x) 





,2,„2/ 



A 

4!< 



d„P" + 2R"[-d lt tp(x) + P„] + m 2 <p{x) + -ip z {x) = 0, 
fyiZ" + 3K»[-dv<p(x) + P^} + 2R^ + X -m 2 + ^ 2 (x) = 0, 
B^K" + AM»[-d^(x) + P M ] + eR^K" + ^ip(x) = 0, 



A 
24 



0, 



(3.46) 
(3.47) 
(3.48) 
(3.49) 
(3.50) 



where terms of power higher than 4 are omitted for convenience. By supposing the self-consistent canonical 
equations 



d»ip(x) = P", 



cq. (3.47) gets satisfied, and eq. (3.46) is transformed to be 

d„S*" = dM^P" - - \™?V 2 (?) - - A = if* (<p{x)) . 

Eqs. (3.48), (3.49) and (3.50) arc simplified to be 



+ 2R fl R» + X -m 2 + ^ 2 (x) = 0, 
d„K^ + QR„K^ + ^<p(x) = 0, 
d^M" + + 8M„iJ" + A = 0. 



(3.51) 
(3.52) 

(3.53) 

(3.54) 
(3.55) 
(3.56) 



The exact solutions of these equations are difficult to derive. However, the behavior of solutions around small 
r is enough for our purpose. For small r, we can replace ip(x) with y>(z) in eqs. (3.54), (3.55) and (3.56). Then 
we can get the following series solutions for P M , and M M 



V 11 + 0(r 3 ), 



p" + 0(r 3 ), 



1 1 1 

2 r ~ 6 £ 

Ci 1 Ci Ci 2 

— rr -\ — oj to r 

r 3 2 r 8 

± -I — i J £ _| i 

2 r 5 r 4 2 r 3 3 r 2 

2 ' \ 1 2lj 2 

9Ci + 7C 2 CW + -^ r C 3 p fl - 



U) 

24 



(3.57) 
(3.58) 
(3.59) 



rp^ + Oir 2 ), 
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where C2 = G\ + , u> = y ^- + and C3 is a constant. While the series expressions of S*^ and 

can be derived from equations (3.51) and (3.53) 

,2 



= ip{z) 



-^ + 0(r 3 ), 



-m 2 - ^ip 2 {z) 



-p" + 0(r 3 ), 



(3.60) 
(3.61) 



where a = p^p^. Here we notice that the foregoing solutions are very similar to that which we obtained in 
section 2.2.1. 



3.2.2 Derivation of the Discretized Lagrangian 

Having obtained the solutions in the last subsection, it is straightforward to derive the discretized Lagrangian 
following the procedures in section 3.1.3. By the definition of eq. (3.26), taking limits term by term as we did 
in section 2.2.2, then we can get the lattice Lagrangian 

(3.62) 



-lattice = lim lim lattice = lim. lim V] 

z l ^x z Lp(z*)—xtt(x z ) z z ^x z <p(z 1 )^k£>(x 1 ) — ' 



3 1 



Pip 



i=0 
3 



-imW)-^V)-A 



i=0 
3 



y — x' 



i=0 



2e 2 p. 



1 



6 

m 2 t \(/> 2 (x l ) 
a 



W) - 4>(x*j\ 



2a 



3 3 

1^ iftPVCa*) - 0(^)] 3 - ^ E ^ W) - 

i=0 



24^ a 

i=0 



In derivations of this Lagrangian, we have supposed the symmetrical lattice spacing y l — x l — e, and in the 
limiting procedures above, we have set the constants of integral C2 and C3 to be zeros, so C\ = — 3^2^ , in 
order that the limits can be well defined. From this lattice Lagrangian, further assumptions on infinitesimal 
lattice spacing e will lead to the lattice Lagrangian of interactions as we expected. 



4 Topologically Massive Gauge Theory 

In this section, we discuss the topologically massive gauge theory introduced in [1, 2]. For an Abelian gauge 
field in three dimensional Minkowski space-time supplemented with the Chern-Simons topological term, its 
Lagrangian is given by 

J2f = ^ V + l^ x F^A x , (4.1) 

where the three dimensional Lorentzian metric is 77 — diag(l, — 1, —1), and e^ uX is the totally antisymmetrical 
tensor with the convention e 012 = 1. In three dimensions, A is a minus constant of mass dimension, and 9 is 
dimensionless. In [1, 2], the authors show that the gauge field excitations are massive. The massive excitations 
can be understood as follows. The Eulcr-Lagrangc equation for (4.1) is 

d^ v - 6 -he vafi F a p = 0. (4.2) 
The dual field strength of F a p is defined as 

*F»= l -e^F a p, (4.3) 
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e 2 A 2 



* p a 



then F a p can also be expressed by its dual field strength 

F" v = e^ a *F a . (4.4) 

Substituting this expression into the equation of motion (4.2), we get an equation of motion for the dual field 
strength 

(^ a d / , + 0A7 ? !/a )*F a = O, (4.5) 

which furthermore implies 

(e px „d x - 0A VvP ) (V^ + 9ArT) *F a = d a dp - VaP (d^ 

= - (d^ + 9 2 A 2 ) *F P = 0. (4.6) 

In the calculation above, we have used the Bianchi identity 

d a *F a = 0. (4.7) 

Eq. (4.6) clearly demonstrates that the dual field *F 13 is massive. For another analysis based on the Dyson- 
Schwinger equation, see [13]. 

4.1 De Donder-Weyl Equation for Gauge Theories 

We start with the Lagrangian (4.1) supplemented with the gauge fixed term 

if = ^F^F^ + 6 -e^ x F^ v A x + ±d u A^A». (4.8) 

The reason for supplementing the gauge fixed term will be clear soon later. Performing the covariant Legendre 
transformation 

c) 9^ 1 1 

= -F'" y + -d v A» + -e* vX A\. (4.9) 



d(dM A £ 2 

Without the gauge fixed term, because F^ is a antisymmctrical tensor, we could get a constraint relation for 
the symmetric part of n llv 



T fiV , VII 



(4.10) 



As is well known in conventional Hamiltonian theories, the appearance of the constraint condition implies the 
singularity of the Legendre transformation, which is ascribed to the gauge invariance of the original theories. 
For the gauge fixed Lagrangian (4.8), the constraint condition (4.10) disappears 1 . Here and hereafter we use 
the gauge £ = A, so the derivative of gauge field can be expressed by its conjugate momentum as follows 

d"A" = Att^ - ^6s» vX A x . (4.11) 



The covariant Hamiltonian is derived as 



JT = VA"^ -2 = ^'"V - ^0e^\^Ax + ^0 2 A^. (4.12) 



1 In this paper, we always add a gauge fixed term to cure the problem of irregularity of Lagrangian. Another approach is 
adding Lagrange multipliers to the Hamiltonian to handle the gauge constraints as that did by Dirac and Bergmann [14-16]. 
More geometrical analysis on the constraints have been carried out in [17]. 



- 16 - 



The canonical Hamiltonian equation is given by 

d ^ = U = - \^^ A \ ( 4 - 13 ) 

which together means that the Euler-Lagrange equation is 

M 0M" - A6e a ^d a A f3 = 0, (4.15) 

which is the gauge fixed version of eq. (4.2). We noticed that a mass term appears in the covariant Hamiltonian 
(4.12), which might be ascribed to the massive property of this gauge theory as we discussed in the beginning 
of section 4. However, the contribution of this mass term is canceled by that of the second topological term 
in (4.12) secretly. So the Euler-Lagrange equation is still of the same massless form as that of eq. (4.2). 

Similar to the discussions in section 3.1.1, we designate the Lagrangian to be the independent integral of 
Hilbcrt 



and designate 



dA v > 

then we derive the De Donder-Weyl equation for the Abclian gauge field 



.17) 



8S fJ '\ 

[A ll ,^ = — \=0. ( 4 - 18 ) 



For the covariant Hamiltonian (4.12), its corresponding De Donder-Weyl equation is 

KdS^dS }L _K A A 
2 dA v dA v 2 dA v A 4 



9/ + ^° ,7 - -^0e^ x ^A x + - r 9 2 AnA' J ' = 0. (4.19) 



We employ the "embedding method" to find solutions for the De Donder-Weyl equation (4.19). At first, 
a solution of the gauge fixed Euler-Lagrange equation (4.15) can be given by 

B v (x) = (-^Cr + lj B v {z) + u)Ce va pB a (z)(x fi - sp) - ^u 2 CB a (z){x a - z a )(x v - z v ), (4.20) 

r = (x» - z^){x^ - u) = A9, 

where z M and C are constant. We have used the initial condition B^{x)\ z = B^(z) to normalize the solution. 
Noticing that the De Donder-Weyl equation (4.19) is quadratic, we can suppose a closed quadratic polynomial 
for Si* 

S»(x) = S*»(x) + P^(x) \A v {x) - B v (x)] + Y^(x) [A a (x) - B a (x)} [A p {x) - B p (x)} . (4.21) 

Obviously, Y tJ - a/3 (x) should be symmetrical about its indices a and j3. Substituting this expression into 
eq. (4.19), we can get a quadratic polynomial of \A^(x) — B fl (x)]. Then the De Donder-Weyl equation (4.19) 
can be satisfied by 

- d»B v {x)P^ + ^P^P^ - ^0e^ x P^B x (x) + ^e 2 B^x)B^(x) = 0, (4.22) 

[A a (x] - B a (x)\ N a = 0, (4.23) 
[A a (x) - B a {x)] [A p (x) - B p {x)] M afi = 0, (4.24) 
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where N a and M a " are defined by 



A 



A 



A 



N a = (0 M P" a - 'jee^n^ + ^9' 2 B a (x)j - 2 \d„B p (x) - AP^ + ^9e^ x B x (x)j Y» a0 , (4.25) 

M afi = d^Y^P + 2A?/ A y' JQI % At , - A9^ x s^ a Y^ x + ^9 2 if?. (4.26) 
As that in the "embedding method" , supposing the relations 

dpBpix) = AP^ - ^9e^ x B x (x), (4.27) 

d^ a = ^Oe^ a n^ - ^9 2 B a (x), (4.28) 

eq. (4.23) gets satisfied. Eqs. (4.27) and (4.28) are self-consistent because they are just the canonical Hamil- 
tonian equations (4.13) and (4.14) and B^(x) solves the Euler-Lagrangc equation (4.15). While eq. (4.22) is 
transformed to be 

d„S^(x) = d»B v {x)P»» - ^P^P^u + ^9e^ x P^B x (x) - ^9 2 B^x)B^(x) = X*(b^ 
The solutions for can be obtained from eq. (4.27) 



P, 



1 „ C 2 2 



<)f) 



Oe^ x B A + —luCB^{z)(x v - z v ) - — uCB v {z){x - z)„ 
29 



- -uj 2 CB a {z)(x - z) a e^x{x - z) x - ^wCB a (z)(x - z) a r]^ 
The solution for 5* M can be obtained from eq. (4.29) 



5*^(x) 



1 



-[C z + C\uj9 9uj A C z r + dr-^ 



cxp 



125 

C 2 + ^9u} 3 C 2 I s*exp s ds 



B a {z)B a (z)(x- zf 



B a (z)(x-z) a (x-zY 



While a solution of eq. (4.24) for Y tia ^(x) can be given by 

1 UJ 



Y^(x) = ^v^cot {^k x {x - z) x ^j n^k", a = k x k x , 



where k 11 is a constant vector. 



(4.29) 



(4.30) 



(4.31) 



4.2 Derivation of the Discretized Lagrangian 

Following the procedures in section 3.1.3, we can derive the lattice Lagrangian from the solution (4.21) for 
the gauge field. Because the solutions (4.29), (4.30) and (4.31) seem complicated, the calculations will be 
laborious. A useful skill to perform these calculations is to employ the additive property of limiting, that is, 
we calculate the limits term by term. Supposing the symmetrical lattice spacing y 1 — x 1 — e and e is small 
enough, we can get 

lattice = lim lim JUice = Km lim V gVL__^(g) (4.33) 

2 



i=0 



y — x' 

2 



1 £ sgn^) v Myi) _ Mxl) i t Aa{yl) _ Aa{xl) i B£ l sgn(l)£l ,x \ Myl) _ Mxi) \ Ax{xi) 



2 — e 

i=0 



L?C- ^-ljC) £sgn(i) AM) - Ai(x*) M**) 



10 
99 



i=0 
2 



j*c - -cjc) [My 1 ) - MJj\ ^V)- 

' i=0 
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The constant coefficients C\ and C2 in eq. (4.30) have been settled to be zeros in order to avoid singularities 
in the limits. For infinitesimal e, A a (y 1 ) — A a (x 1 )— >0, we can also make the replacement 



A a (y i ) - A a {x l ) ^ diA a (x). 



Then the lattice Lagrangian in (4.29) will approximate to the continuous one 

& = ^d»A v d»A v + -e^ x d^A„A x , 



(4.33) 



(4.34) 



which is just the gauge fixed Lagrangian (4.8) with the gauge choice £ = A. 

From the foregoing discussions, it seems that we did not derive anything new from the solutions of the De 
Dondcr-Wcyl equation. However, we caution that the lattice Lagrangian (4.32) is not the unique Lagrangian 
we can derive from the solutions of the De Donder-Weyl equation. As we did in section 2.1.3, we can derive 
another lattice Lagrangian from another solution of the De Donder-Weyl equation. This is partially because 
eq. (4.31) is not the only solution of eq. (4.24). A second solution of eq. (4.24) can be given by 



x) = - _L -^=tan 



k\(x 



?7 a/ V. 



(4.35) 



Together with the solutions in eqs. (4.29) and (4.30), eq. (4.35) completes a new solution for the De Donder- 
Weyl equation (4.19). From this new solution, repeating the procedures as we did in deriving the lattice 
Lagrangian (4.32), we can derive a new lattice Lagrangian 



SA 



w 



lattice 



LjC 2 6A a {x)A a {x) - -wC^sgn(i) My 1 ) - A^x 1 ) A t ( 3 



(4.36) 



6^2 -sgn«e wA [A v (y*) - A„(x 1 )} A x (x l ) - -LoC^My 1 ) - M^)} A a (x i ). 

i=0 i=0 

This Lagrangian is similar to the lattice Lagrangian (4.32). Actually, the Lagrangian (4.36) is a part of the 
lattice Lagrangian (4.32). What causes these differences is that the solution (4.31) contributes to the lattice 
Lagrangian in the approximation of small enough e while the solution (4.35) does not. For infinitesimal e, 
making the replacement (4.33), we then derive the corresponding continuous one 



\ + C^j 6e^ x d„A„(x)A x (x) - LuC 2 6A a (x)A a (x), 



where the corresponding Euler-Lagrangian equation will be 



e^ a d v 



2C 2 



2C+ 1 



wtf" )A a (x) =0. 



(4.37) 



(4.38) 



Obviously, eq. (4.38) is extremely similar to eq. (4.5), which is the equation of motion of the dual field strength 
*F fJ -. Actually, when C = 1 + \/3, 2 c+i = 1> ec l- (4-5) is recovered up to a minus sign, which has no effect on 
the equation of motion of second order. Eq. (4.38) demonstrates the massive property of the gauge field as 
that eq. (4.5) does. 

As in section 2.1.3, we can find a third solution of eq. (4.24) for Y^ a ^ 

Y^{x) = -~ (l + u^- tan(a,y^)) rf^(x - z)», r = (x- z^x - z),. (4.39) 
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Associated with cqs. (4.29) and (4.30), eq. (4.39) constructs a third solution for S^. Following the similar 
limiting procedure above, we can derive the third Lagrangian 



££ = -^d Vb A v d»A v 



2C ) 6e^ x d^A v (x)A x {x) - (1 + 2C)Cw9A a (x)A a (x). 



(4.40) 



This Lagrangian has a mass term, but its kinetic has the minus sign, so this is also a ghost kinetic term. This 
is because the solution (4.39) contributes to the Lagrangian opposite to that of the solution (4.31). 



5 Mass Generating Mechanism for a Fermion Field 
5.1 De Donder-Weyl Theory for a Fermion Field 

We consider a massless fermion field in D dimensional Minkowski space-time with the Lagrangian 



(5.1) 



where the Lorentzian metric if = diag(l, — 1, — 1, • • • ) is mostly negative diagonally, and a^ LV = i [7^, 7"] is 
the spin matrix. The last term of (5.1) is introduced in [12] to make the Legendre transformation regular 2 . 
The non-regularity of the Legendre transformation can also be cured by adding the Wilson's term d^rf" d v i\) 
[18, 19] or some more general term as investigated in [20]. However, an interesting feature of the last term in 
(5.1) is that it is a total divergence 

d^a^d^ = ^ [^ v d^\ , (5.2) 
so it has no effect on the equation of motion. The momentums are defined by the Legendre transformations 



7T 



which means the derivative of fields can be expressed with momentums 

i ( „ 11 



A \ n Ta » 2D^l r 



where D is the dimension of space-time, and t m „ is defined by 



D-2 



1 



The covariant Hamiltonian can be given by 



D - 1 



^7 

? (IV • 



(5.3) 
(5.4) 

(5.5) 
(5.6) 

(5.7) 
(5.8) 



A 



ITT^Tuu'K 1 ' H iblb 



A 



i 1 -u 1 i 1 



2D- 1 



2D- 1 



^7 /i 7T A 



From this Hamiltonian, We can obtain the canonical equations of motion for ip and 7r A 



A 



i 1 



2D - 1 



u 1 D 



(5.9) 
(5.10) 



2 Of course, we can also employ the Dirac-Bergrnann method [14-16] to handle the problem of non-regularity. 
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which together means that the Euler-Lagrangian equation is 

i'fd^ = 0. 



(5.11) 



We notice that although a mass term appears in the Hamiltonian (5.8), the Euler-Lagrangian equation is 
massless. The contribution of the mass term is canceled secretly. For -0 and 7^, the canonical equations of 
motion arc 



d<ff 
dip 



2D-1 
1 _ u ID 



2D — 1 



AD -I 



Supposing the Lagrangian is the independent integral of Hilbert, we have 



Furthermore, by designating 



-dS^ dS^ 



dS^ 



dip ' 
dS^_ 

dip 



we then obtain the De Dondcr-Wcyl equation for a fermion field 



dip 



For the Hamiltonian (5.8), its Dc Donder-Weyl equation is 



dip) 



d u S" + A 



dS" dS" 1 D - : 

i—^ — T,, v — — H WW 

dip M dip 40-l w 



+ A 



2D -I \dip 



( 



7m i> ~ V>7^ 



dS^ \ 



0. 



(5.12) 
(5.13) 

(5.14) 

(5.15) 
(5.16) 

(5.17) 
(5.18) 



We employ the "embedding method" to seek solutions of the De Donder-Weyl equation (5.18). At the 
first step, a solution of the Euler-Lagrangian equation (5.11) is given by 



X (i)=^yr 7 / + lJ X (z)+C 7a (x- Z r^), X = k a k a , r = k a (x~z) c 



(5.19) 



where k a and z a are constant vectors, and C is a constant of mass dimension. Here the initial designation 
x(a;)|z = xi z ) is employed. Because eq. (5.18) is quadratical, we can suppose a closed quadratical polynomial 
for S>* 

S»(x) = S*» + P» [i>{x) - X (x)} + $(x) - x(x)] P» + ${x) - X (x)] R» [i>{x) - x(x)) ■ (5.20) 
Substituting (5.20) into (5.18), we obtain a quadratical polynomials about [ip(x) — x{ x )]- The Dc Donder-Weyl 
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equation (5.18) can then be satisfied by supposing 



A 



i 1 



iP^P v + -^-j (P^x(x) - x(xh,Pn + 



1 D 



4D-1 



X{x)xix) 



[il>{x) - X (x)] : - d^xW 1 + A 

[${x) - x(x)] : fyP" - Wd^xix) + A 



i 1 



iP'VJZ" + -— - (P% - x{x)l^) + 



1 D 



4D-1 



Employing the "embedding method" , we can suppose that 

a M x(a;) = a 



i 1 



-d^ = A 



i 1 „ 1 -D , ' 



2D- 1 



4D-1' 



= o, 

(5.21) 
= 0, 

(5.22) 
= 0, 

(5.23) 
= 0. 

(5.24) 

(5.25) 
(5.26) 



then eq. (5.23) can be satisfied. Eqs. (5.27) and eq. (5.27) are just the Hamiltonian canonical equations (5.9) 
and eq. (5.10). Similarly, we can suppose 



duXix) = A 



-fyP" = A 



i i 

iP a T afl - -yj— jX(xh^ 



i 1 s„ ID, ' 
n + Tn-TXW 



2D- 1 



4D-r 



(5.27) 
(5.28) 



then eq. (5.22) gets satisfied. They are just the Hermitian conjugate versions of eqs. (5.25) and (5.26). Based 
on these assumptions, eq. (5.21) is now transformed to be 



= -A 



+ P"a„x(x) + dyxW = J?*[x(x),x(x) 



Substituting the solution (5.19) into (5.25) and (5.27), we can obtain solutions for P M and P M 



C/D 



P p = 



2^ + a(a^"-^ 



r7 a fc Q - 7 a (x - z) Q 
r7 Q fc Q - 7q(x- z) q 



While the solution for 5*^ is given by 



C 2 



S*" = -5(fl-l)rt"x(^). 



A solution for W of eq. (5.24) can be given by 

ID -I 



AD- 2 



wtanh(wr)A: AI , D ^ 2, 



R» = --— — Ark' 1 , 
4D-1 



D = 2, 



(5.29) 

(5.30) 
(5.31) 

(5.32) 

(5.33) 
(5.34) 
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where lo 2 = ^^A 2 . 

5.2 Derivation of the Lattice Lagrangian for a Fermion field 

Following the limiting procedures, we can derive the lattice Lagrangian for the Fermion fields. Furthermore, 
we can derive the continuous Lagrangian as we did. For simplicity, here we only give the continuous one 



—tp^ a k a k' J 'tp — ipj^ip 
A 



(5.35) 



So we obtain a fermion Lagrangian with a mass term and some modifications. For C = 0, we can get a 
conventional Lagrangian for massless fermions, but for C ^ 0, a mass term and modifications appear. We 
might imagine that we can recover the original Lagrangian (5.1) if we employ other solutions of the De Donder- 
Weyl equation as we did for scalar fields and gauge fields. However, it seems it is very difficult to achieve this 
goal. The reason is that the term d^ipa* 11 ' d v ip is difficult to be derived by the limiting procedure as we suggest. 
As compensations for the loss of this term, a mass term and some modifications appear in the Lagrangian 
(5.35). This situation is quite similar to that in the topologically massive gauge theories. In that case, the 
kinetic term disappears while a mass term appears. 

6 SU(2) Yang-Mills Gauge Theory 

6.1 De Donder-Weyl Equation of the Yang-Mills Gauge Theory 

For the SU(2) Yang-Mills gauge field in four dimensional Minkowski space-time, we consider the Lagrangian 
with gauge-fixed term 

Sg = ^F^F a ^ + hSpAovBrA*, F% = d^A a „ - d v A atl + ge abc A b ^A w , (6.1) 

where e abc are structure constants of SU(2) Lie algebra. We introduce a constant A, which should be —1 for 
the four dimensional Lorentzian metric ij = diag(l, —1, —1, —1). The gauge fixed term is introduced in order 
to make the Legendrc transformation regular. The covariant Legendre transformation is given by 

<9_S? 1 1 

Making the gauge choices £ = A, the derivatives of the field can be expressed by the momentum 

0M£ = ATi afW - ge abc A»A v c . (6.3) 
Then the covariant Hamiltonian is found to be 

je = d^A v ^ v |7r a ^7r^ - ge abc A bil A cv ^ + ^e abc A^A cv e afh A)A v h , (6.4) 

which yields the canonical equations of motion 

dp,A av = = k-K ailv - ge abc A bfi A cu , (6.5) 

-d^ v = = \e dbc A hl Ac V s dfa A^ f + e dba A bll (w dflu - n d ^). (6.6) 

OA av A 1 

Together they mean that the Euler-Lagrangian equation 

d^A av + ge abc (d„A^A cu + 2A^A CU - A^A^) + g 2 e abc e cd f A b ^A» = 0. (6.7) 
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Supposing the Lagrangian to be the independent integral of Hilbcrt, we get 

dx» p p dA av 



^= — = d^+d,A a „ — . (6.8) 



By designating 



afiu 



dA a „ 

we then obtain the De Donder-Weyl equation for SU(2) gauge theory 

dS> i \ 



(6.9) 



+ & ( A au ,n a ^ = — = 0. (6.10) 



For the Hamiltonian (6.4). the De Donder-Weyl equation is 



d * 3 ' + imaai - + if - < 6 -"> 

We employ the "embedding method" to seek the solutions for the De Donder-Weyl equation (6.11). At 
the first step, we need a solution for the Euler-Lagrange equation (6.7), one of which solutions can be given by 

Bg(x) = ^^JacobiCN ( r, ) <p%, a = k^k^, r = u)W{x - z) M (6.12) 



V2g \'V2 
where is a constant tensor, which is restricted by the relations 

- Sab, rtkp = 0, (6.13) 

where 8 ab is the Kronccker tensor. A solution of these restrictions can be given by 

p£ = §y + P a d> + 6 3 a d», 6% = -1, c% = -1, d% = -1, (6.14) 
6% = 0, 6% = 0, c% = 0, = 0, fc% = 0, fc% = 0. 

This solution has been listed in [21]. For more solutions of the SU(2) Yang-Mills equation, see [21, 22]. We 
suppose a series solution for 5 M of the following type 

= S*» + P a » a [A aa (x) - B aa (x)} + [A aa (x) - B aa (x)} [A bf} {x) - B bp (x)} (6.15) 

+ Kbc 1 \A aa {x) - B aa (x)} \A bp (x) - B bp {x)] [A cl (x) - B cl {x)] 

+ M abcd X \A aa {x) - B aa (x)] \A bp (x) - B bp {x)\ [A CJ (x) - B cl {x)] [A dX {x) - B dX {x)] 

+ Kb°f d f XT l A ac,(x) - B aa (x)] [A bp (x) - B bp {x)\ \A cl {x) - B cl (x)} \A dX {x) - B dX (x)} [A fT (x) - B fr (x)} 



Substituting this series into eq. (6.11), we can obtain a scries about \A aa (x) — B aa (x)]. Eq. (6.11) can be 
satisfied if we suppose the coefficients of this scries to be zeros. For the term independent of fields, we have 



^P CHLU P a ^ + ^e ahc B hii B cv e^ h B^Bl 



fyS*" = P afia d l _ t B aa + ge abc B bf ,B cu P a ^ - 
For the term of power 1, we have 

d„P bflf3 ~ ge acb B c ^{P atlf3 - P af3 ») + 



3.16) 
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For the term of power 2, we have 



= [A bp {x) ~ B bp (x)] [A cl (x) ~ B cl (x)) Tg, 



T, 



«f 7 + (-d,B aa + AP a , a + ge abc B b ,B ca ){KZ P1 + *LT + <T) 
£ W + 0(V + R caJ) 9S abc P a ^ - 9 e adc B d ,(R^ + R^) 



(6.18) 
(6.19) 



- ge acd B dlM {Rlf + itff 4 ) + f - ^(e^e^ + e hbf e hdc )B d Bj + 2e hbd e hcf B d „B 

4A L ' 

For the term of power 3, we have 

= [A bp {x) - B bp (x)] [A cl {x) - B cl {x)} \A dX {x) - B dX {x)} Q%£, 



f 11 



Qlt = «7 A + (-d,B aa + AP a , a + ge abc B b ,B ca ){MZ P d t ^ bacd - - 



Mi 



fiflyXa 
bcda 



A U V abc + A bac + A bca )\ U ad^ + n da^ v) 9^ { n a d + H da I 



(6.20) 

) 

(6.21) 



jr. 

4A 



„ F a S d Bt (K^ 1 4 4- K^ lX \ - aF adf R f (K Xylil 4 R' A ' 3l/ ' 7 4 K xPiy \ 

gt of^\±\ abc -r J\ bac a- J\ bca ) gt of v \t\ abc -t J\ bac \ b ca > 

~e hdf B]if x + e hfc e hdb B$if~< + e hbf e hdc B} v f3X + e hbc e hfd B p f r? x . 
For the term of power 4, we have 

= [A bp (x) - B b0 (x)} [A cl (x) - B cl {x)] [A d \(x) - B dX (x)} \A fr (x) - B fr (x)] X^ Xt , 

n 2 

y MXt _ p, hrfJ .pi\T , hbc hdf BX 7 t _ np adf i K XTpy . K -X/3r'y , 7^A/3 7 t\ 

A bcdf — U V IU bcdf 4A £ " " 9 £ K^abc + 1X bac + 1X bca > 

i AD i ^a,bcu R A I AT^PfXr , M fj,f3a~/XT , 

bed/ 



4- C— R 4-AP 4-riF abc Bi, B VAT mq ^ 7At 4- Ar^PaiXr . jy-fiPjaXr 



AT 



p,j3-f\o>T 
abcdf 



N 



(6.22) 
(6.23) 

abcdf 



') 



A 



4_ —(K^Pl _i_ K^ pv "l 4- K^ P1V \<K At 4- ft' A T -I- ft' Ar "l 
+ n l-^a&c +iV bac + 11 bca IK^adf^v + lx daf^v ^^dfan v) 



A V M abcd + iU bacd + M bcad + IU bcda )\ K af^ + K fa^u) 



rfiaPfX , M n0ajX , trupjaX , ^jVPfXa 
l abcd lu bacd ' 111 bead 

C/fc £>h^\m abcd -t m bac d ^ IU bcad ^ IV1 bcda I 



ge^ h B hv {Ml 



bed 



bacd 



bead 



M, 



bcda 



For simplicity, we do not display the terms of higher power here. As that in the "embedding method", we can 
make some assumptions to solve these equations. We suppose that 



■Sj.pW = -ge acb B cll {P a ^ - P af}f *) 



—e £ J B cv B f B d , 



3.24) 
3.25) 



then eq. (6.17) can be satisfied. Eqs. (6.24) and (6.25) are just the canonical equations (6.5) and (6.6). They are 
self-consistent because B aa is a solution of Euler-Lagrangian equation (6.7). On these assumptions, eqs. (6.18), 
(6.20) and (6.22) are simplified largely. While cq. (6.16) is transformed to be 



dfiS *» = P a^g^ Baa _ ±pw Pativ + ge abc B bft B cu P a ^ - £- £ abc B b ,B w e afh B» B£ = J?* (B aa (x)) 



(6.26) 



The solutions for S"* M and P a f lv therefore can be obtained straightforwardly. However, the solution for R^ b 
is not easily found. We can suppose the solution for R^ b ^ of the type 



R^ = X 



(r)5 abV af} kv + R(r)e abf (^r?^ - rfif^ + F(r)<^^ + Z(r)^ a tfk» + K{r)5 aW a } ^ 



3.27) 
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Substituting this expression into (6.18), and employing the relations in eq. (6.13), we can then know that 
(6.18) is satisfied if the following ODEs hold 

uja^ + 2A (X 2 - 2R 2 ) + 4 5 0(r)i? - ^ 2 (r) 

W(J ^ + 2A (Z + 2X - K) R - gMr){Z + X - K) + 

dr 2A dr 

dY 

ioa— + 2A (-2Z + 2X - 2K — 3Y) Y — 2g4>{r)R 
dr 

ua^L _ 4A (R 2 + KZ - XZ) + 
dr y ' 2 A 

ujo— - 2A (K 2 - 2R 2 - 2KX + Z 2 ) + 2g6(r)R - ^r0 2 (r) 
dr 2A 

where 4>{r) = ^^JacobiCN (r, -^=J . The exact solutions of these ODEs are difficult to obtain. However, 
from the discussions in sections 3 and 4, we know that only the behavior of the solutions at small enough r is 
important. So we turn to seek solutions at small enough r. In the foregoing ODEs, we designate 

' ' > ^JacobiCN (r, ±= ) (6.33) 



0, 


(6.28) 


0. 


(6.29) 


0, 


(6.30) 


0, 


(6.31) 


0, 


(6.32) 



y/2g \ ' V2J V^g ' 

that is, we replace 4>{r) with Then we can obtain the solutions at small enough r 

A^coth(^), Z=-^coth(^), (6.34) 

F = -^ 1 5 , R=0, K = ^coth(^), 

3A sinh(2^) +Cisinh 2 (-^) 4A v ^ h 

or we can get another solution 

A^tanh(^), ^-^tanh(^), (6.35) 



Y=-^4- ' 2 . J* = 0, A-^tanh( — ). 

3A sinh(2^) + C 2 cosh 2 (-^) 4A 



These solutions have expressions similar to those that we obtained in sections 2 and 3. The solutions for K^bf 1 
and M^ c ^' yX are more difficult to obtain. We expect that they might not contribute to the lattice Lagrangian 
in the infinitesimal lattice spacing limit as happened in sections 2 and 3 for nonlinear theories. 

6.2 Derivation of the Lattice Lagrangian 

In this section, we plan to derive the lattice Lagrangian for the SU(2) Yang-Mills gauge theory. However, we 
are immediately met with problems. The problem is the components of the solution (6.12) are subjected to 
the constraint condition (6.13). If we take the limits A afl (x)—^B afl (z), then the components of A a)Ji (x) must 
also be subjected to the constraint condition (6.13) in order to avoid singularities in the limits. Therefore, 
what we can derive by the foregoing limiting procedures is a Lagrangian with its field components subjected to 
be the constraint condition (6.13). An approach to overcome this problem is to find more general solutions of 
the Euler-Lagrangian equation (6.7). This means we need find solutions with more general initial conditions. 
The initial value problem of the Yang-Mills equation has been discussed in [23] based on theories of partial 
differential equations, and also can be addressed in the De Donder-Weyl theories as it has been did in [24]. So 
it seems we have no problems in finding solutions of more general initial conditions. However, even if we can 
obtain a general solution for the Euler-Lagrangian equation, it is still a challenge to find a manifestly Lorentz 
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covariant solution of the De Dondc-Wcyl equation, due to the difficulty in solving tensor equations. From our 
discussions on the restricted case, it seems that it is possible that we can obtain several kinds of solutions of 
the De Donde-Weyl equation. 

In section 4, we showed that we can derive a massive Lagrangian of first order for the Abclian topologically 
massive gauge theory. Furthermore, for the non-Abelian Yang-Mills theory, the massive infrared behavior 
of the gluon propagator has been observed in large volume lattice simulations [25-27]. So an interesting 
question is whether we can also derive a massive Lagrangian for SU(2) Yang-Mills theory in our present 
framework. The answer is inconclusive. The derivation of lattice Lagrangians depends on the solutions of the 
De Donder-Weyl equation. We can speculate that there might exist a solution which would lead to a massive 
Lagrangian for SU(2) Yang-Mills, but we did not find such a solution so far because of the complexities of tensor 
equations. The solution we can get is a restricted case as we showed in the last section. To derive a massive 
Lagrangian or to develop some qualitative criteria which judge whether we can derive a massive Lagrangian 
is an interesting question which deserves our future considerations. However, our present framework can still 
help us to understand nonpcrturbative behaviors of Yang- Mills theory from another perspective. In [28-30], 
the authors proposed that the 0(3) nonlinear a model 

vn ^ 

£eff = — (^n) 2 + ~(n • d^n x d u ii) 2 + V(n), (6.36) 

where n is a three dimensional unit vector, could be related to the effective low energy Lagrangian of a SU{2) 
Yang-Mills theory, because the Lagrangian (6.36) supports solutions of closed knotted solitons [31-33]. In [34], 
supposing a special gauge field decomposition 

A M = C^n + c^n x n + pd^n + ad^n x n, (6.37) 

the authors show that an effective Lagrangian similar to (6.36) can be derived from the Yang-Mills action 
after one loop functional integration. In our present framework, we can try to find a solution by supposing 
the gauge field expansion (6.37), then by employing the "embedding method" and the limiting procedures, we 
can speculate that we can also derive a Lagrangian similar to (6.36). Therefore, our present framework may 
provide another approach to understand the infrared behavior of Yang-Mills theories, although it is difficult 
to implement. 

7 Further Discussions 

7.1 Discussions on Solution Dependence 

The main points of our paper can be summarized as follows. Starting with solutions of the Hamilton- Jacobi 
equation or the De Donder-Weyl equation, we can derive a lattice Lagrangian by the limiting procedures 
introduced in sections 2 and 3. The derived lattice Lagrangian can be employed to formulate a lattice definition 
of path integrals. These steps make it possible to derive quantum theories from solutions of the Hamilton- 
Jacobi equation 3 . Here problems arise because we can obtain several kinds of solutions of the Hamilton- Jacobi 
equation (or De Donder-Weyl equation for field theories) generally. We can derive several different lattice 
Lagrangians from these different kinds of solutions. Some Lagrangians have the problem of stability, as we 
obtained ghost Lagrangians in sections 3 and 4. There exist natural criteria for mechanics to select out the 
physical Lagrangian but they do not apply to field theories. 

For the harmonic oscillator, we obtain 3 solutions (2.15), (2.16) and (2.18) for its Hamilton- Jacobi equation. 
The solution (2.15) leads to the conventional Lagrangian by our limiting procedures, while solutions (2.16) 
and (2.18) do not. However, we can use boundary conditions to fix the solutions. For a time evolution from t a 

3 For conventional lattice gauge theories, the lattice Lagrangian employs the gauge link variables so the problem of gauge 
invariance is resolved. While our derived lattice Lagrangian depends on the field variables directly, so the gauge invariance is 
broken and we need introduce ghost fields to cure the problem of renormalization. 
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to tb, we have a orbit from (t a ,x a ) to (tb,%b) in the configuration space. The classical action along this orbit 
can be evaluated, which is given by 

S(x b ,t b ;x a ,t a ) = — - — -[cosuj(t b - t a )(xl + x 2 a ) - 2x b x a ]. (7.1) 

Z SlIlLJ[t5 — t a ) 

By designating the constants of integral appropriately, solutions (2.16) and (2.18) can both lead to classical 
actions equivalent to (7.1). So the uniqueness of the problem of mechanics can be resolved with the help of the 
classical action. However, this method cannot be applied to field theories, obviously. The action of field depends 
on the volume of space-time. There is no unique definition for the action of field between two different space- 
time points. We might expect the integrability conditions introduced in section 2.1.1 to restrict the solutions 
of the De Donder-Weyl equation, but they cannot remove the nonphysical solutions. At least they do not 
work for the free scalar field theory. Another possible criterion is the requirement of employing a "complete 
solution" of the De Donder-Weyl equation to derive the lattice Lagrangian. A "complete solution" should 
depend on the n x d constants of the integral, where n is the number of real components of field variables 
and d is the number of space-time dimensions. But this criterion does not apply to the free scalar field theory 
neither, and there are also difficulties in finding "complete solutions" for other field theories. The solutions 
for the gauge field and fermion field given in previous sections are not "complete solutions" . These problems 
imply our procedures are not complete and we need additional conditions to remove the nonphysical solutions. 
However, as we just discussed, it seems no appropriate criteria emerge naturally. 

7.2 comparisons with conventional lattice gauge theories 

Another problem we should mention is how to implement gauge symmetries in our derivation of lattice La- 
grangians. In conventional lattice gauge theories, the gauge link variables U^(x) are employed; therefore, the 
gauge invariancc is manifest in appropriately chosen lattice Lagrangians. In our present framework, we meet 
the same problems that we usually meet in the framework of canonical quantization: The Lagrangian is not 
regular due to the gauge invariance. In this paper, we solve this problem by adding gauge fixed terms. The final 
Lagrangians we derived are gauge fixed Lagrangians but not gauge invariant ones. This is an unsatisfactory 
aspect of our framework compared to the conventional lattice gauge theories. It is a challenge to figure out 
whether we can implement gauge invariancc by appropriate modifications in our present approach. Leaving the 
gauge symmetry problem aside, it seems that our approach has advantages on dealing with infrared behaviors 
of field theories. In section 4, we derived a massive Lagrangian for a topologically massive gauge theories. In 
section 5, we suggested a mass generating mechanism for fermions. In section 6, we discussed the possibilities 
of deriving an effective Lagrangian for a Yang-Mills theory. These advantages are achieved by employing 
special solutions of Dc Donder-Weyl equations cleverly. However, we should caution that these advantages are 
not ensured, because we can always derive physical Lagrangians associated with nonphysical ones while we do 
not have any natural criteria to select out physical ones. These aforementioned potential advantages could be 
verified until we can achieve a better understanding about the problem of solution dependence as we discussed 
in the last subsection. 

8 Conclusions 

The construction of covariant Hamilton- Jacobi theories or De Donder-Weyl theories can be traced back to the 
efforts to formulate the Lorentz covariant canonical quantization from the Dirac canonical bracket 4 . These 
efforts met problems partially because the covariant Hamiltonian docs not correspond to the total energy 
of the system. Rather than attempting to construct a covariant canonical formulation, we suggest deriving 
lattice Lagrangians from solutions of De Dondcr-Wcyl equations by appropriate limiting procedures. It turns 
out we can obtain several different sectors from different kinds of solutions generally. We can obtain a massive 

4 For several suggestions recently, we refer to [35-38] . The De Donder-Weyl equation for gravity of Einstein-Hilbert action has 
been derived in [39]. The application of Dc Donder-Weyl theories to canonical quantum gravity has been suggested in [40, 41]. 
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Lagrangian of first order for a topologically massive gauge theory. We also find a surface term that can 
produce masses for fermions, which is a mass generating mechanism similar to that of topologically massive 
gauge theories. The problems are that we also obtain ghost Lagrangians, which means we might need an 
mechanism to get rid of these nonphysical results. For nonlinear theories like Yang-Mills theories, our analysis 
on restricted solutions suggest we can derive nontrivial lattice Lagrangians as we did for topologically massive 
gauge theories, but these remain speculative until we can handle complicated tensor equations. 
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